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A model  is  developed  wliich  yields  the  best  system  design  (number 
and  location  of  repair  depots,  number  of  repair  channels  at  each  depot, 
and  size  of  the  spares  inventory  at  each  depot)  to  service  a widely  dis- 
persed population  of  repairable  units,  subiect  to  random  failure.  Con- 
sidered in  the  design  model  is  a service  level  constraint  in  terms  of 
spares  availability,  as  well  as  purchase  costs  of  spares  and  repair  chan- 
nels, fixed  costs  of  operating  repair  depots,  and  unit  repair  and  trans- 
portation costs.  Failure  and  repair  times  are  assumed  to  be  exponential 
and  finite  source  queueing  theory  coupled  with  a branch-und-bound  algo- 
rithm are  utilized  in  developing  the  model. 
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1 . Introduction 

A population  of  units  (aircraft  engine  components,  computer  modules, 
lathes,  etc.)  that  randomly  fail  but  are  completely  repairable  requires 
both  spares  and  repair  facilities  to  provide  for  certain  desired  levels  of 
service.  For  example,  it  may  be  desired  to  have  the  system  perform  so  that 
the  probability  that  at  least  B percent  of  the  population  operating  at 
any  given  time  is  1-a  . Another  measure  of  system  performance  often  used 
Is  the  percentage  of  requests  for  spares  met  Immediately  from  on-hand  in- 
ventory. This  latter  criterion  is  referred  to  as  spares  availability  and 
is  the  criterion  used  in  this  model;  that  is,  the  complement  of  the  percent- 
age of  requests  back  ordered. 

The  physical  system  under  consideration  is  of  a cyclical  nature  in 
that  there  is  a fixed  number  of  components,  say  N . If  we  desire  M of 
those  to  be  operating  at  a given  time,  then  there  are  N - M = y spares  in 
the  system.'  Any  component  can  be  in  one  of  three  states:  (1)  operating, 


*We  assume  a component  to  he  "basic,"  that  is,  upon  failure,  the 
entire  component  (not  parts  of  it)  is  replaced  by  an  identical  spare  unit. 
No  indenture  is  considered  in  this  study,  as  opposed  to  the  MOD-METRIC 
model  (see  Muckstadt  (1973)1. 
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(2)  spare,  or  (3)  in  the  repair  facility  (either  being  worked  on  or  In 
a repair  queue). 


In  order  to  determine  the  appropriate  number  of  spares  (we  assume 
M is  given  as  input  and  it  is  desired  to  find  y ) and  the  repair  facility 
capacity  necessary  to  achieve  a specified  service  level,  a stochastic  model 
is  required  since  the  failure  and  repair  processes  are  random.  The  real 
stochastic  process  is  extremely  complex;  many  simplifying  assumptions  must 
be  made  to  obtain  analytical  results.  Therefore,  the  stochastic  process 
modeled  here  is  an  approximation  to  the  true  underlying  stochastic  process, 
and  the  simplifying  assumptions  are  duly  noted.  We  believe,  however,  that 
the  approximation  is  an  adequate  one  and  that  the  simplifying  assumptions 
are  justified,  as  discussed  below. 

2.  Mathematical  Model 

Letting  y equal  the  number  of  spares  in  the  system  and  c equal 
ttie  number  of  simultaneous  repair  channels  ir  the  repair  depot,  we  desire 
to  find  the  values  for  c and  y such  that 

Minimize  ':(c,y)  = C^(c,y)  + Cjj(c,y)  (1) 

c,y 

Subject  to  2 q 1-ot  , (2) 

n=0  " 


where 


C^(c,y)  = annually  prorated  one-time  costs  such  as  purchase 

costs,  salvage  values,  set-up  costs,  etc. 

Cjj(c,y)  = annually  recurring  costs  such  as  repairmen  salaries, 
holding  costs  of  spares,  unit  repair  costs,  etc. 

q = probability  that  n units  are  in  the  repair  queue 

” or  in  repair  when  a request  for  a spare  occurs  (an 

item  fails). 

The  problem  then  is  one  of  finding  the  combination  of  c and  y that 
minimizes  annual  equivalent  costs  over  the  life  of  the  system  subject  to 
providing  a service  level  of  l-CK  (percentage  of  requests  for  spares  filled 
immediately  from  on-hand  inventory). 
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y-1 

Treating  Equation  (2)  first,  7 q is  the  probability  that  when 

-n 

n=i) 

a request  occurs,  at  least  one  spare  is  available;  that  is,  fewer  than  y 

units  are  in  or  awaiting  repair.  Considering  the  repair  process,  q are 

n 

simply  the  conditional  (upon  a failure)  probabilities  that  the  system 
state,  which  is  measured  by  the  number  in  the  repair  system  either  in  queue 
or  being  worked  on,  is  n . These  q^  are  functions  of  the  general  time 

system-state  probabilities  (call  them  p ) which,  with  some  assumptions, 

n 

can  be  determined  from  classical  queueing  theory.  This  will  be  elaborated 
on  below;  however,  first  we  treat  Equation  (1)  In  detail. 


Suppose  the  system  is  expected  to  be  used  for  k years,  and  that 
r is  the  annual  Interest  charge  on  capital.  Then  a one-time  expenditure 
at  the  beginning  of  the  first  year  can  be  converted  to  an  annual  equivalent 

k k 

expenditure  by  applying  the  capital  recovery  factor,  r(l+T)  /[(1+r)  -1] 
[see,  for  example,  Maynard  (1956),  pp.  7-66  to  7-77].  Salvage  value  re- 
ceived at  the  end  of  k years  can  be  converted  to  an  annual  equivalent 

return  (negative  cost)  by  applying  the  sinking  fund  factor,  r/[(l+r)  -1]  . 
Using  the  following  additional  notation. 


'p,y 


P,c 


C 


s,y 


CR 


SK 

R(c,y) 


purchase  cost  of  a spare  ($  per  unit) 

purchase  cost  of  a repair  channel  (S  per  channel) 

salvage  value  of  a spare  after  k years  ($  per 
unit 

salvage  value  of  a repair  channel  after  k years 
($  per  channel) 

operating  cost  of  a repair  channel  ($  per  year  per 
channe 1 ) 

carrying  and  handling  costs  of  a spare  ($  per  year 
per  spare  in  system) 

unit  repair  cost  ($  per  unit  repaired) 
capital  recovery  factor 
sinking  fund  factor 

average  number  of  units  repaired  per  year, 
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Equation  (1)  expands  to 


Minimize  Z(c,y)  = (Cp  - C F^p  + Cp)y 

c,y 

^ ^S.c'^CR  - S.c^SF  + 


(3) 


+ CpR(c,y)  , 


where  R depends  on  c and  y and  the  stochastic  process  model. 


3.  Stochastic  Process  Model 


The  stochastic  process  modeled  here  is  a finite  population  of  items 
that  randomly  fail.  Upon  failure,  the  failed  item  is  removed,  transported 
to  the  repair  facility,  repaired,  and  dispatched  to  the  spares  pool.  Si- 
multaneously a spare  is  requested  and  "plugged  in"  Immediately  if  one  is 
available.  If  no  spares  are  available  ( y or  more  units  in  or  awaiting 
repair)  the  request  is  back  ordered.  The  service  level  constraint  requires 
the  percentage  of  requests  back  ordered  to  be  £ a • This  process  is  sche- 
matically shown  in  Figure  1. 

Assuming  that  failures  and  removal,  transportation,  and  repair  times 
are  random  yields  a cyclic  queueing  process,  since  we  also  assume  complete 
repair  capability  and  no  addition  or  withdrawal  of  items.  When  failures 
are  Poisson  and  removal,  transportation,  and  repair  times  exponential, 
analytical  results  can  be  obtained  [see  Gross  and  Ince  (1977)1.  A further 
simplifying  assumption  in  which  removals  and  transportation  times  are  insig- 
nificant or  are  Included  in  the  repair  times  so  that  the  cyclic  queue  con- 
sists of  only  two  stages — the  ready  unit  stage  and  the  repair  stage — reduces 
the  model  to  that  of  the  classical  machine  repair  problem  with  spares  [see 
Gross,  Kahn,  and  Marsh  (1977)]. 

The  key  assumptions  necessary,  enabling  use  of  the  classical  queueing 
theory  to  determine  q^  and  R in  Equations  (2)  and  (3),  are  then 

1.  Poisson  failure  process, 

2.  exponential  repair  (transport  and  removal)  times, 

3.  all  units  Identical, 
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4.  continuous  operation  of  units  and  repair  facilities, 

5.  steady  state. 

We  now  discuss  these  in  detail. 

3.1  Poisson  failures 

A Poisson  failure  process,  or  equivalently,  exponential  times  to 
failure,  is  commonly  assumed  when  modeling  lifetime  distributions  of  many 
components,  and  especially  those  with  considerable  electronic  hardware. 

It  has  also  been  used  extensively  when  dealing  with  aircraft  engines,  par- 
ticularly jet  aircraft.  While  there  is  great  impetus  to  use  this  assump- 
tion for  mathematical  convenience,  its  final  justification  must  rely  on 
how  well  it  approximates  the  real  world.  In  many  cases,  data  have  borne 
out  the  reasonableness  of  the  Poisson  assumption;  however,  data  are  not 
always  available,  for  example,  for  a new  component.  In  such  cases,  either 
data  on  similar  components  (justifying  the  assumption  for  gas  turbine  ship 
engines  using  aircraft  jet  engine  data,  for  example)  or  characteristics  of 
the  Poisson/exponential  process  itself,  such  as  the  Markovian  property,  the 
long  tall  of  the  exponential  distribution,  and  so  on,  must  be  utilized. 

If  one  considers  the  shape  of  the  exponential  density  function, 
the  long  tall  to  the  right  Indicates  that  one  would  occasionally  experience 
a very  long  time  to  failure,  but  that  most  of  the  failure  tlvies  are  of  mod- 
erate length.  The  Markovian  or  constant  failure  (hazard)  rate  property 
says  that  the  item  exhibits  no  wear-out;  that  is,  an  item  that  has  been 
operating  a long  time  has  no  more  chance  of  failing  than  a brand  new  item. 
Tliese  two  characteristics  seem  to  be  quite  reasonable  for  the  kinds  of  items 
encountered  in  a repairable  Inventory  system — jet  engines,  avionics  equipment, 
etc . 


Another  property  of  the  exponential  distribution  is  that  the  minimum 
of  a set  of  exponential  random  variables  is  also  exponential.  Thus,  if  a 
unit  is  rather  complex  with  many  "pieces,"  each  of  which  can  fail  approxi- 
mately exponentially  causing  the  unit  to  fail,  the  unit  has  an  exponential 
failure  distribution. 
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3.2  Exponential  repair  (tranjport  and 
removal)  times 

In  order  to  use  classical  queueing  theory,  repair  times  must  be 
exponentially  distributed  unless  the  repair  capacity  has  "ample  servers;" 
that  is,  no  queueing  ever  occurs  [see  Posner  and  Bernholtz  (1968)].  The 
same  holds  true  for  remo”-'.'!.  and  transport  times.  It  is  rare  for  repair 
depots  to  have  an  "ample"  number  of  channels;  that  is,  a potential  channel 
for  every  item  in  the  system.  The  ample  server  assumption,  however,  may 
be  more  realistic  for  the  transport  and  removal  functions. 

Unless  there  are  ample  servers,  the  justification  of  exponent iallty 
for  these  times  must  be  done  in  a manner  similar  to  that  used  to  justify 
failure  times.  If  there  is  a sizable  human  element  in  the  service  process, 
one  might  expect  more  symmetrical  distributions  for  service  times  with  in- 
creasing hazard  rates;  i.e.,  the  longer  an  item  has  been  in  service,  the 
greater  the  probability  of  its  being  completed  in  the  next  Instant  of  time. 

However,  if  service  is  mainly  diagnostic  rather  than  routine,  or  if  it  is 
very  heterogeneous  in  nature  (each  repair  differing  from  the  last),  then 
the  exponential  assumption  with  its  long  tail  and  Markovian  property  could 
be  quite  realistic.  Again,  either  a data  analysis  or  knowledge  of  how  the 
process  is  likely  to  work  would  be  required  for  final  justification. 

3.3  Identical  units,  continuous  operating, 
and  steady  state 

These  three  assumptions  are  treated  in  some  detail  in  Barzlly, 

Cross,  and  Kahn  (1977),  and  are  shown  to  be  adequate  for  repairable  item 
systems.  We  discuss  those  briefly  here  and  refer  the  interested  reader  to 
the  reference. 

] 

To  use  classical  cyclic  or  machine  repair  queueing  theory,  all 
units  must  have  identical  failure  and  repair  cliaracterist ics , so  that  it 
is  necessary  only  to  keep  track  of  how  many  are  at  each  stage  in  the  cycle, 
and  not  which  Individiial  units  are  where.  In  reality,  tills  is  generally 
not  so;  some  repaired  items  may  not  be  as  good  as  original  items  and  some 
may  be  better  due  to  technological  learning.  If  items  are  introduced  over 
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time,  later  Items  may  be  more  reliable.  Newer  Items  may  also  be  easier 
to  repair.  To  account  for  this  while  using  classical  queueing  theory, 
we  averaged  characteristics  of  all  items  in  the  population  and  assumed  a 
population  of  identical  items,  all  operating  with  these  "average"  charac- 
teristics. An  analysis  in  the  reference  cited  above  showed  that  unless 
item  characteristics  differed  by  an  order  of  magnitude,  the  averaging 
assumpt'on  was  adequate. 

Classical  queueing  theory  assumes  that  all  ready  items  in  the  popu- 
lation (except  spares)  operate  continuously,  as  do  the  repair  f t ' ' llties. 
Often  in  the  real  situation,  while  both  repair  facilities  and  operating 
units  keep  the  same  hours,  all  population  units  may  not  be  operating  simul- 
taneously. For  example,  a fleet  of  aircraft  with  M total  engines  requires 
a support  system  (repair  facilities  and  spares),  but  not  all  aircraft  operate 
all  the  time.  Furthermore,  those  not  operating  at  any  time  cannot  be  con- 
sidered to  be  spares,  for  they  may  be  regularly  or  randomly  scheduled  for 
operation.  This  is  usually  accounted  for  by  "stretching  out"  the  mean 
time  to  failure  of  each  unit  and  "pretending"  all  units  operate  at  a re- 
duced failure  rate.  Another  adjustment  method  suggested  in  Barzily  £t  al . 
(1977)  is  to  find  an  average  percentage  of  the  population  which  is  operating 
at  any  time  and  to  adjust  M by  this  percentage  while  using  the  true 
failure  rate.  These  two  methods  were  compared  and  found  to  be  in  very 
close  agreement. 


The  problem  of  steady-state  results  was  also  analyzed  in  Barzily  ^ 
al . , where  it  was  found  tiiat  the  steady  state  was  quickly  approached  in 
most  cases  unless  population  characteristics  were  changing  very  rapidly. 
When  designing  support  facilities  for  a population  of  units,  interest  is 
centered  on  the  population  once  it  has  reached  maturity,  rather  than  during 
its  infant  years  when  things  are  changing  rapidly.  Tliat  is,  facilities  are 
designed  to  handle  the  estimated  steady-state  situation,  and  while  the  es- 
timates of  some  steady-state  input  parameters,  such  as  final  attainable 
reliability,  may  be  less  certain,  these  can  be  handled  by  a type  of  sensi- 
tivity analysis  that  considers  different  possible  steady-state  situations. 
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Once  the  above  assumptions  are  made  and  their  ramifications  kept 
in  mind,  a classical  queueing  model  can  be  used  to  yield  the  and 

R of  Equations  (2)  and  (3),  respectively.  We  refer  the  reader  to  Gross, 
Kahn,  and  Marsh  (1977)  and  to  Gross  and  Ince  (1977)  for  details.  The  best 
combinations  of  c and  y that  minimize  the  objective  function  Z(c,y) 
of  (3)  can  be  found  using  a two-variable  search  or  can  be  closely  esti- 
mated using  a heuristic  algorithm  presented  in  Gross,  Kahn,  and  Marsh 
(1977).  ine  model  is  further  extended  below  to  consider  multiple  repair 
depot  locations  as  well  as  simultaneously  provisioning  for  more  than  one 
type  of  item,  which  can  share  repair  depots  but  requires  dedicated  repair 
channels  within  the  depot. 

4 . Multiunit,  Multilocation  Design  Model 

The  basic  model  (s),  described  above  for  a single  population  and 
single  repair  facility,  can  be  coupled  with  a previously  developed  multi- 
activity,  multifacility  design  model  that  utilizes  a branch  and  bound  algo- 
rithm [see  Plnkus,  Gross,  and  Sol  and  (1973)]  to  determine  the  number  of 
spares  and  the  number  of  repair  channels  for  each  type  of  unit  at  each 
repair  facility,  as  well  as  how  many  and  where  the  repair  facilities  should 
be  located  among  available  sites.  We  first  illustrate  the  concepts  of  the 
model  on  a small  problem  in  which  complete  enumeration  can  be  used  in  place 
of  the  branch  and  bound  algorithm,  and  then  illustrate  the  power  of  the 
branch  and  bound  methodology  on  a large  problem. 

Consider  a hypothetical  fleet  of  200  domestic  aircraft  with  four 
types  of  repairable  units  that  can  be  accommodated  in  the  same  repair  depot 
as  long  as  separate  repair  channels  are  used  for  each  of  the  four  unit 
types.  Each  aircraft  requires  three  each  of  unit  types  1 and  2 and  two 
each  of  unit  types  3 and  A.  For  example,  one  unit  type  might  be  a major 
engine  component,  another  a landing  gear  assembly,  and  so  on.  Since  on  the 
average  only  half  the  planes  operate  at  any  one  time,  our  effective  fleet 
size  is  100.  Letting  M^  denote  the  operating  population  size  of  unit 

type  1,  we  have  M^  = M^  = 300  and  M^  = M^  = 200  . We  have  three  possible 
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repair  depot  designs  to  consider.  Design  1 consists  of  one  central  repair 
depot  in  the  midwest  that  supports  the  entire  fleet.  Design  2 consists  of 
two  additional  sites  on,  say,  the  east  and  west  coasts,  so  that  each  depot 
serves  one-third  of  the  fleet.  Design  3 has  two  sites  in  addition  to  that, 
or  five  depot  sites  in  all,  which  are  distributed  around  the  country  so  that 
each  services  one-fifth  of  the  fleet.  No  transshipping  is  allowed;  that 
is,  each  depot  services  a specific  portion  of  the  fleet  and  it  is  assumed 
that  each  site  has  the  potential  to  support  a depot  of  any  desired  capacity. 
Figure  2 shows  the  possible  designs  schematically. 

Given  any  design,  it  is  known  what  portion  of  the  fleet  is  served 
by  each  depot,  so  that  the  previous  model  [the  model  giving  the  solution 
to  Equations  (2)  and  (3)]  can  be  used  separately  for  each  unit  type  popu- 
lation at  each  depot  for  each  design.  We  use  the  acronym  SASPRO  (Spare  and 
Server  Provisioning)  to  refer  to  this  model.  For  example,  consider  Design 
2 (three  depot  locations)  for  unit  type  3.  Each  of  the  three  depots  serves 

a population  of  size  67  ( ^ [100  x 2]  ).  With  the  appropriate  costs, 

failure  rate,  and  turnaround  (removal  + transport  + repair)  time,*  the 
number  of  repair  channels  and  spares  for  each  depot  can  be  determined 
from  SASPRO,  along  with  total  costs,  from  Equation  (3).  To  compare  differ- 
ent designs  for  a given  unit  type,  we  must  account  for  a reduction  in  repair 
time  due  to  closer  depots,  and  must  also  Introduce  transportation  costs, 
which  were  ignored  in  the  basic  model  since  it  had  always  been  assumed  that 
all  units  went  to  a central  depot.  Further,  the  fixed  cost  of  a depot  must 
be  considered  when  the  possibility  of  various  design  configurations  exists. 

Basically,  the  design  model  considers  four  additional  factors  that 
must  be  traded  off:  (i)  the  decrease  in  transportation  costs  when  more 
depots  are  involved,  (ii)  the  reduction  in  turnaround  time  due  to  the 
proximity  of  depots,  (ill)  the  increase  in  fixed  facility  costs  when  depots 
are  added,  and  (Iv)  the  possible  economy  of  scale  (fewer  total  system  spares 


I 


i 

i 


’We  use  in  these  examples  the  machine  repair  (two-stage  cyclic  queue) 
model,  so  that  the  service  operation  is  considered  to  be  one  operation  that 
includes  removal,  transport,  and  repair. 
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and  repair  channels  may  result  from  greater  centralization  of  depots; 
that  is,  a few  large  depots  rather  than  many  small  ones). 

First,  we  consider  only  the  costs  shown  in  Equation  (3)  for  all 
possible  designs  for  each  unit  type.  We  assume  a 20-year  system  life  and 
an  annual  interest  charge  of  10%.  The  desired  spares  availability  (l-a) 
is  90%.  The  failure  rate,  A , is  In  units  of  failures  per  day,  while 
the  mean  turnaround  time  is  in  units  of  days.  All  costs  are  in  thousands 
of  dollars.  Table  I shows  the  results  of  applying  SASPRO  to  all  possible 
designs  for  each  type  of  unit.  Note  that  the  mean  repair  (turnaround) 
time  decreases  as  more  depots  are  added  because  of  the  smaller  transport 
distance.  It  was  assumed  that  transport  time  was  20%  of  the  total  turn- 
around time,  and  it  was  cut  proportionately  when  additional  depots  were 
introduced.  The  costs  shown  in  the  next  to  last  column  are  only  the  costs 
given  by  (3)  that  were  incurred  at  each  depot.  Thus,  for  a design  with 
three  depots  the  costs  must  be  tripled,  as  must  c and  y when  compared 
with  the  single  depot  case.  The  costs  for  the  total  population  are  shown 
in  the  last  column.  The  economy  of  scale  is  rather  evident  here  when  the 
total  system  spares  and  repair  channels  required  for  a single  depot  are 
compared  to  those  required  for  the  multidepot  case.  The  last  column 
clearly  reflects  this. 

If  transportation  time  were  a significantly  large  portion  of 
total  turnaround  time,  the  economy  of  scale  of  a centralized  depot  could 
be  lost.  However,  even  without  this,  when  transportation  costs  are  con- 
sidered, a multilocation  depot  design  could  be  preferable.  Table  11 
introduces  transportation  costs. 

While  transport  time  was  cut  in  proportion  to  the  number  of  depots 
(that  is,  n depots  cut  the  transport  time  portion  of  the  turnaround  time 
to  1/n  ),  the  costs  would  not  be  reduced  by  the  same  magnitude.  We  assume 
here  that  for  three  depots,  transportation  costs  are  cut  in  half,  and  for 
five  depots,  the  transit  costs  are  cut  by  two-thirds.  Furthermore,  we 
assume  that  the  transport  cost  is  a certain  percentage  of  the  purchase 
cost  of  an  item,  but  that  this  percentage  varies  depending  on  the  particular 
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SPARES,  REPAIR  CAPACITY,  AND  COSTS  (10  DOLLARS)  PER  DEPOT  AREA 
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unit  type.  The  percentages  used  are  shown  In  Table  II.  Note,  for  example, 
that  unit  type  1 may  be  transported  quite  inexpensively  with  respect  to  Its 
value  (purchase  cost);  namely,  its  transport  cost  is  only  2%,  while  the 
transport  cost  of  unit  type  2 is  relatively  expensive,  or  23%. 

The  expected  number  repaired  at  each  depot  per  year,  R , comes  out 
of  the  SASPRO  model.  This  figure,  multiplied  by  transportation  costs  and 
number  of  depots,  yields  the  total  transportation  cost  per  year  per  popu- 
lation of  each  type  of  unit  for  each  possible  repair  depot  design.  Adding 
this  to  the  costs  given  in  Table  I yields  what  we  refer  to  as  the  total 
variable  cost  associated  with  each  :'nit  type  population  and  depot  design, 
and  is  given  in  the  last  column  in  Table  II. 

If  the  fixed  costs  of  operating  a repair  depot  are  ignored,  by 
searching  the  final  column  of  Table  II  we  see  that  the  best  design  for 
each  unit  type  is  as  follows: 

Unit  Design  (Number  of 
Type  Number  depots) 


11  1 

2 3 5 

3 2 3 

4 2 3 

When  introducing  the  fixed  costs  of  operating  a depot,  however,  it  may  be  ad- 
visable to  select  a design  other  than  that  which  appears  best  when  only  the 
variable  costs  are  considered.  Table  TIT  displays  both  variable  and  fixed 
costs.  The  circled  elements  indicate  the  best  designs  when  only  the  vari- 
able costs  are  considered.  Table  IV,  however,  shows  by  complete  enumeration 
the  total  costs;  that  is,  the  variable  plus  the  fixed  costs  for  all  designs. 
Note  that  Design  2 Includes  Design  1 and  that  Design  3 Includes  Designs  1 
and  2.  In  other  words,  if  we  decide  upon  Design  2,  then  if  the  variable 
costs  of  a particular  unit  type  are  less  for  Design  1,  we  would  use  only 
one  depot  for  that  particular  component.  Thus,  when  considering  which  is 
the  best  design,  when  examining  Design  2 we  can  choose  the  lowest  variable 
cost  between  Designs  1 and  2,  and  when  considering  Design  3,  we  can  select 
the  lowest  variable  cost  of  all  three.  Table  IV  indicates  the  best  solution 
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to  be  Design  2 (three  depots)  where  unit  types  2,  3,  and  4 are  to  use 
three  depots  (one-third  of  the  population  serviced  by  each),  and  unit 
type  1 to  use  only  the  central  depot;  that  is,  the  entire  population  is 
to  be  sent  to  one  central  depot.  Here  the  economy  of  scale  for  number 
of  spares  and  channels  outweighs  the  transportation  savings. 

For  larger  problems,  complete  enumeration  is  not  possible,  but  the 
branch  and  bound  algorithm  referred  to  previously  and  described  in  Pinkus 
et  al . (1973)  is  quite  efficient.  We  present  an  illustration  in  the  fol- 
lowing section. 

4.1  Use  of  the  Pinkus  ££  aj[.  branch  and 
bound  algorithm 

In  general,  if  there  are  m allowable  designs,  the  total  number 

m 

of  cases  that  must  be  tried  to  completely  enumerate  a solution  is 

x=l 

2™-!  , since  all  combinations  of  designs  taken  one  at  a time,  two  at  a time, 
...,  m at  a time  must  be  looked  at.  If  there  are  n possible  locations, 

the  maximum  number  of  allowable  designs  is  2^^-!  (any  location  can  be  in 
or  out  of  tlie  design)  so  that  the  total  number  of  cases  to  consider  for 
complete  enumeration  is 


^ ev)  = - > 


In  the  example  to  follow,  there  are  seven  allowable  designs  which  would 
require  2^-1  = 127  cases  to  evaluate. 


In  the  preceding  example,  which  was  solved  by  complete  enumeration, 

3 

there  were  three  possible  designs  yielding  a total  of  2-1=7  cases. 
Although  it  appears  in  Table  IV  that  only  three  cases  are  being  evaluated. 
Column  2 is  really  two  cases  (Design  2 and  Design  2 or  1)  and  Column  3 is 
really  four  cases  (Design  3,  Design  3 or  1,  Design  3 2,  and  Design  3 

or  2 or  1 ) . In  certain  cases,  when  one  design  is  embedded  within  another. 


18  - 


T-367 





I 

1 

hL 


the  larger  always  dominates,  making  it  unnecessary  to  actually  evaluate 
the  smaller.  For  example,  since  Design  1 is  embedded  within  Design  2, 
Design  2 or  1 would  always  be  less  expensive  than  Design  2 by  itself  since 
the  fixed  costs  are  the  same  for  both,  but  the  variable  costs  might  be 
less  as  some  unit  types  might  prefer  Design  1 while  others  might  prefer 

Design  2.  Thus,  in  reality  there  may  be  fewer  than  2"'-!  cases  to  con- 
sider, and  when  all  combinations  of  designs  involving  n locations  are 
allowable  there  are  considerably  fewer  cases  required  to  be  searched  than 
the  upper  bound  given  by  (4).  Because  of  the  embedding  of  designs,  it 
turns  out  that  the  number  of  cases  to  be  searched  is  never  more  than 

2 -1  . Hence  a better  upper  bound  to  the  number  of  cases  required  to  be 

searched  when  there  are  m (£  2''-l)  allowable  designs  is  mln{2"'-l,  2'^-l} 

The  branch  and  bound  algorithm  structure,  under  the  worst  case, 

would  require  searching  2’^^^-!  possible  nodes  [see  Plnkus  £t  (1973)  , 
for  details].  However,  the  algorithm  in  all  problems  tried  so  far  has 
needed  to  search  only  a small  portion  of  these  prior  to  finding  the  optimal 
design.  The  use  of  this  algorithm  is  illustrated  by  the  following  example. 

Consider  a support  system  for  repairing  and  storing  the  spares  of 
16  unit  types.  Seven  different  designs  for  the  system  have  been  proposed, 
with  at  most  nine  depots  being  considered  for  the  entire  system.  These 
seven  designs  range  from  a completely  centralized  system  (Design  1) , which 
Includes  only  one  depot,  to  a decentralized  system  (Design  7)  with  nine 
depots.  The  depots  Included  in  each  design  are  given  in  Table  V. 


TABLE  V 

SYSTEM  DESIGNS  FOR  PROBLEM  CONSISTING 
OF  16  UNIT  TYPES 


Design  Number 

Depot  Numbers 

1 

5 

2 

3,7 

3 

1,5,9 

A 

2, 4, 6, 8 

5 

2, 4, 5, 6, 8 

6 

1,3, 4, 6, 7,9 

7 

i 

1,2, 3, 4, 5, 6, 7, 8, 9 
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The  16  unit  types  to  be  repaired  under  this  system  have  varying 
population  sizes,  failure  rates,  average  repair  times,  and  costs  associated 
with  their  repair,  the  purchase  and  storage  of  spares,  the  purchase  of  repair 
channels,  and  travel  to  depots  for  repair.  Although  too  numerous  to  in- 
clude here  in  their  entirety,  these  parameters,  along  with  the  required 
number  of  servers  and  repair  channels  given  by  SASPRO,  are  presented  in 
Table  VI  for  each  unit  type  if  repaired  under  Design  1.  The  entire  set 
of  parameters  were  considered  for  all  seven  designs  in  turn,  for  all 
unit  types,  in  order  to  determine  which  of  the  seven  designs  minimized 
the  total  variable  cost  of  the  repair  support  system  for  each  unit  type. 

The  results,  that  is,  the  "best"  design  for  each  unit  type,  are  given  in 
Table  VII.  However,  using  these  designs  requires  a system  that  consists 
of  all  nine  depots,  resulting  in  an  additional  cost  (fixed  cost)  of  $49 
million  (see  Table  VIII).  This  brings  the  total  cost  of  the  system,  that 
is,  the  sum  of  the  fixed  and  variable  costs,  to  $246.6  million. 

Applying  the  branch  and  bound  algorithm  to  minimize  the  sum  of  the 
variable  and  fixed  costs  of  this  repair  support  system  produces  the  system 
design  and  variable  costs  given  in  Table  IX.  Vith  the  exception  of  unit 
types  1,  4,  and  9,  the  variable  costs  are  higher  than  those  shown  in  Table 
VII.  Thus,  13  of  the  16  unit  types  are  not  being  repaired  under  the  de- 
signs that  minimize  their  individual  variable  unit  costs.  However,  the 
overall  system  design  found  by  the  branch  and  bound  algorithm  has  reduced 
the  fixed  cost  of  operating  depots  to  $23.5  million  (only  five  of  the  nine 
depots  are  required),  resulting  in  a total  system  cost  of  $242.6  million — 
a saving  of  $4  million. 

Thirty-three  iterations  of  the  branch  and  bound  algorithm  were  re- 
quired to  solve  this  problem.  This  took  11.52  seconds  on  an  IBM  370/148 
computer. 

4.2  Further  work 

Two  areas  of  further  work  are  being  considered.  First,  we  would 
like  to  incorporate  space  constraints  into  the  branch  and  bound  algorithm. 
These  constraints  would  enable  limitations  to  be  placed  on  the  total  amount 
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TABLE  IX 

BRANCH  AND  BOUND  SOLUTION 


Unit  Type 

Design  Number 

Variable  Cost 
($  X 1,000) 

1 

1 

15,747.2 

2 

4 

4,743.2 

3 

5 

21,214.3 

4 

1 

8,466.5 

5 

5 

19,389.6 

6 

1 

6,865.1 

7 

4 

16,808.7 

8 

4 

13,758.8 

9 

1 

14,585.3 

10 

5 

11,358.4 

11 

5 

16,033.4 

12 

1 

12,346.3 

13 

5 

16,473.6 

14 

4 

7,823.1 

15 

4 

21,225.5 

16 

5 

12,287.1 

Total  Variable  Cost: 

$219,126.1 

of  space  available  at  each  depot  for  repair  channels  and  the  storage  of^ 
spares.  At  present  the  design  model  assumes  that  a depot  can  handle,  from 
a space  point  of  view,  anywhere  from  one  to  k unit  types,  where  k in 
the  above  problem  was  16. 

Second,  we  are  investigating  the  possibility  of  performing  a sensi- 
tivity analysis  of  designs  for  unit  types,  in  terms  of  their  various  pos- 
sible parameters.  Hopefully,  this  would  lead  to  a narrowing  of  the  system 
designs  that  would  need  to  be  considered  for  large-scale  problems. 
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